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Stability analysis of collective neutrino oscillations in the supernova accretion phase 

with realistic energy and angle distributions 
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We revisit our previous results on the matter suppression of self-induced neutrino flavor con- 
versions during a supernova (SN) accretion phase, performing a linearized stability analysis of the 
neutrino equations of motion, in the presence of realistic SN density profiles. In our previous numer- 
ical study, we used a simplified model based on an isotropic neutrino emission with a single typical 
energy. Here, we take into account realistic neutrino energy and angle distributions. We find that 
multi-energy effects have a sub-leading impact in the flavor stability of the SN neutrino fluxes with 
respect to our previous single-energy results. Conversely, realistic forward-peaked neutrino angular 
distributions would enhance the matter suppression of the self-induced oscillations with respect to 
an isotropic neutrino emission. As a result, in our models for iron-core SNe, collective flavor conver- 
sions have a negligible impact on the characterization of the observable neutrino signal during the 
accretion phase. Instead, for a low-mass O-Ne-Mg core SN model, with lower matter density profile 
and less forward-peaked angular distributions, collective conversions are possible also at early times. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Supernova (SN) neutrinos are important astrophysical 
messengers to probe the flavor mixing in unique condi- 
tions [l|,[2|. I n particular, renewed attention is being paid 
to collective features of flavor transformations j3Hl5| in- 
duced by v-u self-interactions [16h18| in the deepest stel- 
lar regions, near the neutrino-sphere. The observed col- 
lective phenomena of synchronized [19( and pendular |20| 
oscillations, and the splits [2l| - |23j in the observable SN 
neutrino spectra are subjects of intense investigations 
(see [24J for a recent review). 

The implicit assumption in the characterization of the 
self-induced oscillations is related to the flavor evolution 
in the deepest SN regions, being driven by the only large 
neutrino densities n v . However, during the SN accre- 
tion phase (at post-bounce times i p b < 0.5 s) also the net 
electron density n e is expected to be large. As pointed 
out in |25j , when n e is not negligible with respect to 
n v , the large phase dispersion induced by the matter for 
z/s traveling in different directions, would partially or 
totally suppress the collective oscillations through pecu- 
liar "multi-angle effects". At this regard, in two previ- 
ous papers [26|, |27J we performed a numerical study of 
the matter suppression during the accretion phase, char- 
acterizing the matter and the neutrino density profiles 
with the results from recent long-term SN hydrodynami- 
cal simulations [28[ , based on three flavor Boltzmann neu- 
trino transport in spherical symmetry. Multi-angle nu- 
merical simulations of the non-linear neutrino equations 
of motion in the presence of large neutrino and matter 
densities are computationally very demanding. Indeed, 
these equations are "stiff" and their solutions generally 
involve a fast-changing combination of multi-frequency 
oscillations. In order to reduce this numerical complex- 



ity, we assumed a simplified neutrino emission model. 
Namely, we considered all the neutrinos to be emitted 
"half-isotropically" (i.e. with all outward-moving angu- 
lar modes equally occupied and all the backward-moving 
modes empty) from a common neutrino-sphere with a 
single representative energy In this framework, we found 
that the presence of a dominant matter term inhibits the 
development of collective flavor conversions. The mat- 
ter suppression ranges from complete (when n e 3> n v ) 
to partial (when n e >n u ) suggesting, in principle, time- 
dependent features. In particular, for the iron-core SN 
models we analyzed, we found complete matter suppres- 
sion of the collective oscillations for post-bounce times 
t p b<0.2 s and i p b>0.4 s and partial flavor conversions 
(with electron antineutrino survival probability P ee ~ 
0.5) in the intermediate time range. 

Our works have stimulated further independent inves- 
tigations of these effects [29M31J . In particular, in [30| it 
has been proposed to study the matter suppression of the 
self-induced oscillations applying the linearized stability 
analysis of the neutrino equations of motion, recently 
worked out in |32| . This method is particularly useful to 
circumvent the challenges of full numerical multi-angle 
simulations for the flavor evolution. These are not nec- 
essary if one is interested only in the issue of the flavor 
stability of the dense neutrino gas. Indeed, the stability 
analysis would allow one to determine the possible on- 
set of the flavor conversions, seeking for an exponentially 
growing solution of the eigenvalue problem, associated 
with the linearized equations of motion for the neutrino 
ensemble. With this approach, taking as benchmark for 
the SN density profiles and for the neutrino emission a 
15.0 M Q iron-core SN model from the Garching group, 
the stability analysis has shown a complete suppression 
of the collective oscillations for all the duration of the ac- 
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FIG. 1: Supernova neutrino flux spectra (in arbitrary units) 
for 10.8 Mq progenitor at t p b = 0.225 s (left panels) and 
for 8.8 Mq progenitor at £ p b = 0.250 s (right panels). The 
energy (upper panels) and the angular (lower panels) spectra 
are shown for v e (dotted curves), 77 e (continuous curves) and 
v x (dashed curves) (see the text for details). 



cretion phase [301 ] . Subsequently, the same technique has 
been applied to our numerical results |31| . finding per- 
fect agreement with our sequence of complete and partial 
matter suppression. 

Motivated by these interesting papers, we find useful to 
improve our previous results applying the stability analy- 
sis to our models, including also realistic energy and angle 
neutrino distributions. The plan of our work is as follows. 
In Sec. 2 we introduce the SN models used to character- 
ize the SN densities and, the neutrino energy and angle 
fluxes. In Sec. 3 we introduce the setup for the flavor- 
stability analysis, describing the non-linear equations for 
the neutrino flavor evolution in SNe, and the consistency 
equations coming from their linearization. In Sec. 4 we 
present our results for the stability analysis of the matter 
suppression during the accretion phase for two different 
SN progenitor masses. Finally, in Sec. 5 we comment on 
our results and we conclude. 



II. SUPERNOVA MODELS 

We considered the core-collapse supernova simulations 
of massive progenitor stars with 8.8 and 10.8 M© progen- 
itor from Ref. [28j , taken as benchmark for our numerical 
study in [27|. The first SN belongs to the class of O-Ne- 
Mg-core progenitors [33|, |34[ and represents the thresh- 
old between thermonuclear explosions and core-collapse 
supernovae [28l . |35| . The second is an iron-core progen- 
itor [3g. All models were evolved consistently through 



core collapse, bounce and the early post-bounce phase 
up to several seconds after the onset of explosion j28| . 
We remind that in order to trigger the explosion for 
the 10.8 M Q progenitor model, the heating rates have 
been artificially enhanced in the gain region where neu- 
trinos deposit energy in order to revive the stalled bounce 
shock. 

In our previous study, wc schematically took a mono- 
energetic v ensemble, with a representative energy E = 
15 MeV. Neutrinos of different species were assumed 
to be emitted half-isotropically by a common spherical 
"neutrino-sphere," in analogy with a blackbody emis- 
sion. However, realistic supernova simulations show 
that v angular distributions at the decoupling are far 
from being half-isotropic and are flavor-dependent (see, 
e.g., |30j,l37(). In order to fix a common neutrino-sphere 
radius r = R for the flavor evolution, consistently with 
our choice in [27| we take the radius at which the v e 's 
angular distribution has no longer significant backward 
flux, i.e. a few % of the total one. This typically is in 
the range R ~ 50 - 100 km (see Fig. 4 in (27|) ■ It is 
convenient to parameterize every angular mode in terms 
of its emission angle Or relative to the radial direction 
of the neutrino-sphere. For a half-isotropic distribution 
the occupation numbers are distributed as dn/d cos Or = 
const, or cquivalcntly the radial fluxes are distributed as 
dF/d cos Or oc cos Or [38[ . A further simplification is ob- 
tained if one labels the different angular modes in terms 
of the variable u = sin 2 Or, as in [32|, [38|. Note that for 
a half-isotropic emission at the neutrino-sphere the v an- 
gular distribution of the radial fluxes is a box spectrum 
in < u < 1, since dcos0R/du oc (cos^) -1 cancels the 
cos Or dependence previously mentioned. 

From SN simulations of |28| we extract the angle and 
energy distributions F Ua (E, u) of the different neutrino 
species. In Fig. 1 we show the (angle-integrated) flux 
energy spectra F Va (E) (upper panels) and the (energy- 
integrated) flux angular spectra F Va (u) (lower panels) for 
u e , v e , v x , where this latter indicates the non-electron fla- 
vors. Here, we are showing these fluxes for the 10.8 M Q 
SN progenitor at t p b = 225 ms (left panels), and for the 
8.8 M SN progenitor at t p b = 250 ms (right panels). 
The angular variable < u < 1 is based on R = 69 km 
for the 10.8 M Q model, and on R = 47 km for the 8.8 M 
model. The energy and angular distributions are normal- 
ized to the total neutrino number fluxes of the different 
species (in arbitrary units in the Figure). 

Concerning the neutrino energy distributions for the 
different flavors, the angle-integrated e nergy spectra are 
well-represented by the renowned form 39|, |40j 
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where for a given flavor L v is the neutrino luminosity, 
(E v ) is the average energy and j3 is the energy-shape 



parameter, defined as |39|, |4 



III. SETUP OF THE STABILITY ANALYSIS 



A. Equations of motion 



P = 



(e 2 ) {E v y- 



(2) 



i.e. it is a dimcnsionlcss parameter containing informa- 
tion on the second moment of the distribution, (E 2 ). In 
general, L U: (E v ) and (E 2 ) are all functions of time, and 
are extracted directly from the simulations at hand (see 
Fig. 2 in [II]). In the case of 10.8 M Q we find a typical to- 
tal number flux hierarchy expected during the accretion 
phase, numerically F Ve : F Ve : F Vi = 1.20 : 1.00 : 0.34. 
The first part of the hierarchy is caused by the dclcp- 
tonization of the collapsed core whereas the second is 
caused by the absence of charged current interactions for 
neutrino species other than v e and v e . Passing now to 
the case of the low-mass 8.8 M Q SN progenitor, we realize 
that the spectrum of the non-electron flavors is less sup- 
pressed with respect to the electron (anti)ncutrino spec- 
tra, namely F Ue : F Ve : F Vx = 1.23 : 1.00 : 0.72. This dif- 
ferent behavior is due to the fact that in low-mass core 
SNe, the duration of the accretion phase is very short 
(i-e. i p b <0.03 s in Fig. 2 of Ref. [111). At later times, 
the neutrino luminosities can be approximated by diffu- 
sion and the spectra are determined from thermalization 
processes (e.g., inelastic scattering on electrons/positrons 
and elastic scattering on nucleons). Since these neutral- 
current processes do not distinguish the different flavors, 
this leads to reduced spectral differences between elec- 
tron and non-electron species (see, e.g., |4l|). 

Considering now to the angular flux spectra, we re- 
alize that in the case of the 10.8 M Q model, these are 
significantly forward enhanced (i.e. peaked at small u) 1 
with respect to the half-isotropic emission model. Re- 
markably, the angular distributions of the three v species 
at the conventional neutrino-sphere are rather different: 
since the v x 's decouple at smaller radii with respect to 
i/ e 's and F e 's, their distributions are more suppressed in 
the direction tangential to the neutrino-sphere (i.e. at 
u ~ 1). In the case of 8.8 M Q the angular spectra for 
the electron species are less forward-enhanced than in 
the previous case, and with less pronounced differences 
with the non-electron species. As for the energy spectra, 
this behavior is due to the dominant role of the neutral- 
current processes in the spectra formation. We also note 
that for the models we considered the angular distribu- 
tions do not cross at any u. 1 



1 As recently discussed in 1421 . possible crossings in the angular 
distributions of different flavors could lead to an enhancement of 
the flavor instability. 



We work in a two-flavor oscillation scenario, associ- 
ated to the atmospheric mass-square difference Am 2 tm = 
2 x 10 -3 eV 2 and and with the small (matter suppressed) 
in-medium mixing Q c g = 10 -3 [43|]. We will always as- 
sume inverted mass hierarchy (Am 2 tm < 0) where v's can 
exhibit flavor instabilities for the flux ordering present 
during the accretion phase [6[. For this flux hierarchy, 
three-flavor effects, associated with the solar mass split- 
ting are negligible 0. Following [32[, we write the equa- 
tions of motion for the flux matrices §e,u as function of 
the radial coordinate. The diagonal ^>e,u elements are 
the ordinary number fluxes F Va {E, u) integrated over a 
sphere of radius r. We normalize the flux matrices to the 
total V e number flux n Pc at the neutrino-sphere. Conven- 
tionally, we use negative E and negative number fluxes 
for anti-neutrinos. The off-diagonal elements, which are 
initially zero, carry a phase information due to flavor 
mixing. Then, the equations of motion read |32|, Ejj 



id r $ E ,u = [H E ,u,®E,u] 



(3) 



with the Hamiltonian \li 



mm 



1 / M 2 \ 

H E ,u = — (■— + y/2G F Ni 
v u \2E J 

y/2G F 



Airr 2 



+oo pi 

dE' / du 

-oo JO 



1 - v u v v 



$ 



E'.u' 



•(4) 



The matrix M 2 of neutrino mass-squares causes vacuum 
flavor oscillations. The matrix Ni = diag(n e ,0,0) in fla- 
vor basis, contains the net electron density and is re- 
sponsible for the Mikheyev-Smirnov-Wolfcnstcin (MSW) 
matter effect [Iff with the ordinary background. Finally, 
the term at second line represents the v-v refractive term. 
In particular, the factor proportional to the neutrino ve- 
locity v u = (1 — uR 2 /r 2 ) 1 / 2 (38[ in the v-v interaction 
term implies "multi-angle" effects for neutrinos moving 
on different trajectories [5|, [l~8|, |4J: Ea] ■ In order to prop- 
erly simulate numerically this effect one needs to follow 
a large number [O(10 3 )] of interacting neutrino modes. 



B. Stability conditions 

In order to perform the stability analysis we closely fol- 
low the prescriptions presented in [321 ] and summarized 
in the following. At first we switch to the frequency vari- 
able u) = Aml tm /2E so that E(u) = \Aml tm /2u\ and 
we introduce the neutrino flux difference distributions 



g u _ u = g(u,u) defined as 



|Am 2 



9w,u 



2u 2 



x\e(w)[F Ve (E(w),u)-F Vm (E(w),u)] 



(5) 



+e(-u)[F Vm (E(u),u)-Fv B (E(w),u)] 



normalized to the total v e flux at the neutrino-sphere. 
Then, we write the flux matrices in the form [32f 
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where Tr$ w . u is conserved and then irrelevant for the fla- 
vor conversions, and the initial conditions for the "swap- 
ping matrix" in the second term on the right-hand side 
are s UjU = 1 and S Ut u — 0. Self-induced flavor transi- 
tions start when the off-diagonal term S uv , grows expo- 
nentially. 

In the small-amplitude limit |>S^,ii| <C 1, and at far 
distances from the neutrino-sphere r 3> R, the linearized 
evolution equations for S UtU in inverted mass hierarchy 
(Am 2 tm < 0) assume the form [32| 

id r Su,u = [w + u(A + e/i)]S' w , u 

- fj, I du'(ko'(u + u')g u >, u >Su>,u> , (7) 

where 



du du> g u 
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quantifies the "asymmetry" of the neutrino spectrum, 
normalized to the total v e number flux. The v-v interac- 
tion strength is given by 
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while ordinary matter background term is given by 



A = V2G 



1?- 



F n e 



2r 2 

10 ! (Y±- 

f!5 



10 10 g/ cm3 / V 10 km 



R 



where Y e is the net electron fraction, and p is the mat- 
ter density. The radial distance r is expressed in km, 
while the numerical values of /i and A in the two previ- 
ous equations are quoted in km -1 , as appropriate for the 
SN case. 

One can write the solution of the linear differential 
equation [Eq. J?])] in the form S u ,u — Qu,v,e~ with 
complex frequency fl = 7 + in and eigenvector Q UtU . A 
solution with k > would indicate an exponential in- 
creasing S Ut u, i-e. an instability. The solution of Eq. Q 
can then be recast in the form of an eigenvalue equation 
for Qu^u- Splitting this equation into its real and imagi- 
nary parts one arrives at two real equations that have to 
be satisfied 321 
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A flavor instability is present whenever Eqs. © admit a 
solution (7, k). 



IV. APPLICATION TO OUR SUPERNOVA 
MODELS 

In this Section we will present our results for the stabil- 
ity analysis of the self-induced flavor conversions for two 
representative SN simulations based on different progen- 
itor masses. 



A. 10.8 M@ progenitor mass 

We start our investigation with the case of the 10.8 M Q 
iron-core supernova. We briefly remind the reader our 
previous results presented in Ref. [27|. For this SN model, 
the net electron density n c and the neutrino densities n„ 
for different post-bounce times were shown in Fig. 5 of 

The numerical results of the 



27 



our previous paper 
multi-angle flavor evolution for our schematic model with 
neutrinos of different species all emitted half-isotropically 
with a single energy, were shown in Fig. 7. In particu- 
lar, it was represented the survival probability of electron 
antincutrinos P ee , at different post-bounce times in the 
presence of matter effects and for n e = 0. The matter 
effects produced a complete suppression of the flavor con- 
versions for t p b ,< 0.2 s and i p b > 0.4 s (where n e /n Pe » 1, 
see Fig. 8 in 27J) and partial flavor conversions (with 
electron-survival probability P ee ~ 0.5) at the intermedi- 
ate times (where n e /np c > 1). 

We now compare these numerical results with what 
we find using the stability analysis. For the same post- 
bounce times considered in our previous paper, we show 
in Fig. 2 the radial evolution of the eigenvalue n deter- 
mined from the solution of Eqs. (J9)). We consider the 
following cases: (a) n e = and a half-isotropic neutrino 
emission (dashed curves) , (b) dense matter effects and a 
half-isotropic neutrino emission (continuous curves) and, 
(c) dense matter effects and non-trivial neutrino angular 
distributions (dotted curves). In all these cases we have 
always assumed quasi-thermal neutrino energy distribu- 
tions, parametrized as in Eq. ([TJ. 

We start discussing our results for the case of n e = 0. 
We realize that when the neutrino system enters into an 
unstable regime (k > 0), the k function rapidly grows 
from zero to a peak value greater than one. Compar- 
ing this result with the numerical solution presented in 
Fig. 7 of [23| , we find perfect agreement between the on- 
set of the self-induced flavor conversions determined nu- 
merically and the position of the peak in the k function. 
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FIG. 2: 10.8 M© progenitor mass. Radial evolution of the 
K function at different post-bounce times with n e = for a 
half-isotropic neutrino emission (dashed curves) and in pres- 
ence of matter effects, with a half-isotropic neutrino emission 
(continuous curves) and with flavor-dependent angular distri- 
butions (dotted curves). Quasi-thermal v energy spectra are 
assumed (see the text for details). 



Remarkably, the multi-energy effects included in the sta- 
bility analysis do not seem to produce major changes in 
the onset of the flavor conversions with respect to what 
observed in our previous mono-energetic investigation. 
Of course, once conversions are triggered, multi-energy 
effects would play a crucial role in determining the split- 
ting features in the final neutrino spectra. For their char- 
acterization, multi-energy and multi-angle simulations of 
the flavor evolution are mandatory (see, e.g., [HI]). 

We now discuss the case with realistic matter density 
profiles and a half-isotropic neutrino emission. As ex- 
pected, the flavor instability is strongly suppressed with 
respect to the previous case with n e = 0. In particu- 
lar, the k function is non-zero only at intermediate post- 
bounce times, i.e. t p h = 0.225,0.3,0.325 s in Fig. 2, con- 
sistently with the occurrence of the flavor conversions 
in [231 ■ I n these cases the rise of the k function is shifted 
at larger radii (by ~ 400 — 500 km) with respect to the 
case with n e = 0. In particular, the instability occurs 
when n e ~ n Pa . Moreover, the peak value of k in these 
cases reaches at most ~ 0.5, implying a slower growth 
of the instability with respect to the case with n e = 0. 
In the intermediate time snapshots, where the flavor con- 
versions are not completely suppressed, we find an agree- 
ment between the onset of Fig. 7 in [27( and the position 



of the peak in the k function. 

Finally, we consider the case in which also the flavor- 
dependent forward-peaked neutrino angular distributions 
are taken into account. We find that the k function 
is further suppressed with respect to the half-isotropic 
case. This is consistent with the expectation that the v-v 
strength is weaker for forward-peaked distributions, mak- 
ing the system more stable under the effect of the matter. 
We also note that at £ p b = 0.325 s when the accretion is 
at the end, the effect of the angular-distributions is less 
pronounced, since these become less forward-enhanced, 
as discussed in Sec. 2. 

In order to validate these findings we also performed 
multi-angle and multi-energy simulations of the flavor 
evolution, for the cases with realistic matter densities 
and neutrino energy and angle distributions. We used 
N u = 1200 angular modes and Ne = 80 energy modes. 
Our main goal is to determine if flavor conversions are 
present or not, without being interested in the exact 
final outcome of the numerical simulations. Therefore, 
we did not require the perfect numerical convergence of 
our results. In the cases in which we have not found a 
positive k with the stability analysis, also the numeri- 
cal evolution has shown no flavor conversion. Moreover 
at tpb = 0.225 s, where k < 0.1, the numerical simula- 
tion gives a complete suppression of the flavor conver- 
sions. Conversely, for the two other post-bounce times 
{tpb = 0.3, 0.325 s) where k reaches a peak between 0.2 
and 0.4, we find partial flavor conversions. However, in 
runs with different numbers of angular and energy modes 
the final value of P ee ranges between ~ 0.8 — 0.9, suggest- 
ing that these partial flavor conversions would be practi- 
cally negligible in the characterization of the supernova 
neutrino signal during the accretion phase. 

Finally, we mention that we applied the stability anal- 
ysis also to the case of a 18.0 M iron-core supernova, 
considered in [22]] • However, since the results obtained 
are similar to the ones discussed in this Section, for the 
sake of the brevity, we do not show them here. 



B. 8.8 Mq progenitor mass 

We now analyze the case of a SN with a low-mass 8.8 
M Q O-Ne-Mg core. We remind that in this case there is 
an absence of an extended accretion phase, since the ex- 
plosion succeeds very shortly after the core-bounce. For 
this model, the net electron density n e and the neutrino 
densities n v at different post-bounce times (i p b < 0.25 s) 
are shown in Fig. 11 of [27]]. The numerical results of the 
multi-angle flavor evolution for our schematic model with 
neutrinos emitted half-isotropically with a single energy 
were shown in Fig. 13 of [271 ] . From this Figure one can 
realize that the matter suppression of self-induced flavor 
conversions is never complete, since the matter density is 
very low with respect to iron-core progenitor. In Fig. 3 
we plot the corresponding radial evolution of the k func- 
tion, for the (a), (b), (c) cases introduced in the previous 



K 



3 
2.5 

2 
1.5 

1 

0.5 





3 
2.5 

K 2 
1.5 

1 

0.5 




Lpb 



0.08 s^ 




Lpb 



0.14 s 



lpb 



0.16sh ~t 



t pb = 0.25 sh 



250 500 750 1000 

r (km) 



250 500 750 1000 

r (km) 



FIG. 3: 8.8 Mq progenitor mass. Radial evolution of the 
K function at different post-bounce times with n e = for a 
half-isotropic neutrino emission (dashed curves) and in pres- 
ence of matter effects, with a half-isotropic neutrino emission 
(continuous curves) and with flavor-dependent angular distri- 
butions (dotted curves). Quasi-thermal v energy spectra are 
assumed (see the text for details). 



Section. In the case of a half-isotropic neutrino emission 
with n e = (dashed curves) we find once more a perfect 
agreement between the peak of the k function and the 
onset of the flavor transitions found numerically. 

We now discuss the matter effects. At first, we real- 
ize that the presence of non-trivial angular distributions 
(dotted curves) does not lead to a further suppression 
of the instability with respect to the case with a half- 
isotropic emission (continuous curves). This different be- 
havior with respect to what we have seen in the case of 
10.8 M© SN model, is consistent with what is shown in 
Fig. 1, i.e. the angular spectra of different flavors for the 
8.8 M Q SN are significantly less forward-peaked than in 
the case of the 10.8 M Q SN. Therefore, their effect on 
the flavor stability is less pronounced. Looking in details 
at the different time snapshots in Fig. 3, we realize that 
at i p b = 0.08 s (when n Pe <n e ) the k function presents 
a peculiar shape coming from the interplay between the 
self-induced and the matter effects with a comparable 
strength. The k curve broadens at smaller r with respect 
to the case with n e — 0. The peak of k is slightly reduced 
by the matter effects. However, its position remains the 
same, giving the onset of the flavor conversions found nu- 
merically. At later times, when n e dominates over n Pc , 
the matter suppression becomes more relevant. However, 
since at most n e > 2n Pc , the suppression is never as strong 
as in the case of the iron-core SNe. Also in this case, the 
matter effect shifts at larger radii the onset of the flavor 
conversions. 



V. CONCLUSIONS 

We have performed a linear stability analysis of the 
self-induced flavor conversions during the accretion phase 
for two SN models with different progenitor masses. We 
characterize the SN densities, the neutrino energy and 
angular spectra, with results from recent SN hydrody- 
namical simulations. We compared this method with 
the numerical results of the flavor evolution presented in 
our previous paper [27] , where we assumed a schematic 
v emission, with that all the neutrinos streaming half- 
isotropically from the neutrino-sphere with a single repre- 
sentative energy. For the case of an iron-core 10.8 M Q SN 
(and 18.0 M Q , not shown) we found that the continuous 
v energy spectra do not play a crucial role in the issue of 
the stability of the neutrino ensemble with respect to our 
previous single-energy results. Conversely, the presence 
of forward-peaked v angular distributions significantly 
reduces the strength of the v-v interaction term. This 
effect would enhance the suppression of the self-induced 
oscillations, due to the dense matter term. We find that 
the flavor instability is completely suppressed for large 
part of the duration of the accretion phase, excepct for 
a small time window around i p b ~ 0.3 s. However, we 
checked with multi-energy and multi-angle simulations 
that the effect of these partial flavor conversions would 
be negligible in the characterization of the observable su- 
pernova neutrino flux. Our result is in agreement with 
the stability analysis recently performed with a 15.0 M Q 
SN model from the Garching group J3fl|. In that paper, 
a complete matter suppression of the self-induced flavor 
conversions was found for all the duration of the accretion 
phase. In the case of a low-mass O-Nc-Mg SN with 8.8 
Mq progenitor, where the accretion phase is extremely 
short, the matter density profile is lower and the v an- 
gular distributions less forward-peaked than in iron-core 
models. As a consequence, we found that also with re- 
alistic angular distributions flavor conversions would be 
possible at early times, i.e. at tpb^O.2 s. The different 
pattern of the self-induced flavor conversions for these 
two SN progenitors could be, at least in principle, a tool 
to distinguish between iron-core and O-Nc-Mg core SNe. 

The matter suppression of the self-induced flavor con- 
versions for iron-core supernovae would imply that the 
original neutrino spectra will be processed only by the or- 
dinary Mikheyev-Smirnov-Wolfenstein effect in the outer 
stellar layers. This effect would allow to distinguish the 
neutrino mass hierarchy through the Earth matter ef- 
fects j47| or the rise time of the SN neutrino signal [48| , 
in the case of a #13 neutrino mixing angle as "large" as 
currently measured by the Daya Bay |4{j and Reno [50| 
reactor experiments. These recent measurement confirm 
and greatly strengthen the significance of earl y hi nts sug- 
gested by the long-baseline v^-v,, experiments [5l|,[52J and 
Double Chooz reactor experiment [5j|, especially when 
analysed in combination with other oscillation data 54 j . 
This intriguing consequence implies that the characteri- 
zation of the v flavor evolution during the accretion phase 



is crucial for the interpretation of the observable SN v sig- 
nal. At this regard in |29J], due to a different choice of the 
neutrino angular distributions, a less pronounced matter 
suppression has been found for iron-core SNe. In future 
it will be mandatory to apply the stability analysis also 
to other SN models in order to understand how generic 
arc the results obtained till now. 
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